Entanglement and non-Markovianity of quantum evolutions 
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We address the problem of quantifying the non-Markovian character of quantum time-evolutions 
of general systems in contact with an environment. We introduce two different measures of non- 
Markovianity that exploit the specific traits of quantum correlations and are suitable for opposite 
experimental contexts. When complete tomographic knowledge about the evolution is available, our 
measure provides a necessary and sufficient condition to quantify strictly the non-Markovianity. In 
the opposite case, when no information whatsoever is available, we propose a sufficient condition for 
non-Markovianity. Remarkably, no optimization procedure underlies our derivation, which greatly 
enhances the practical relevance of the proposed criteria. 



The actual dynamics of any real open quantum system 
is expected to deviate to some extent from the idealized 
Markovian evolution that arises from the conditions of 
weak (or singular) coupling to a memoryless reservoir 
U ■ While quantum optics provides realizations that 
are extremely well approximated by such an evolution, 
soft or condensed matter systems evolve subject to condi- 
tions that are generally unsuited to be treated within the 
Born-Markov framework This is particularly the case 
when considering interacting many-body systems, where 
the subsystem's coupling strength may be comparable to 
the coupling to the bath [7j . The exact details of what 
makes a given quantum evolution non-Markovian may be 
complicated, and in many cases, especially when think- 
ing about many-body systems, an accurate microscopic 
model of the system-bath interaction may actually be un- 
feasible. It would therefore be very useful to define some 
simple measure that captures, in some form, the fact 
that the evolution departs from strict Markovianity. This 
problem was addressed by [B| in the context of abstract 
quantum channels and, very recently, an optimization- 
based measure of non-Markovianity founded upon the 
behaviour of the trace distance under complete positive 
(CP) trace-preserving maps have been proposed in 0. 





FIG. 1: Schematic illustration of the envisaged scenario to 
measure the non-Markovian character of an unspecified dy- 
namical map. An arbitrary quantum system, possibly multi- 
level and with an internal dynamics, is subject to the action 
of a local bath, depicted as a golden glow. The system is 
initially prepared in a maximally entangled state |$) with an 
ancilla which is kept shielded from the bath. The green line 
represents some typical decay of the initial entanglement for 
a Markovian evolution while the red line corresponds to a 
possible non-Markovian decay. In this case, the local action 
of the bath is no longer represented by a continuous family 
of CP propagators © and the entanglement between system 
and ancilla is no longer constrained to decrease monotoni- 
cally. It is this deviation 2> B ' that allows us to estimate the 
non-Markovianity of the process despite ignoring any details 
about the evolution itself. 



Formally, the dynamics of a quantum system given by 
a family of trace-preserving maps £ T is called Markovian 
if it defines a one-parameter semigroup of CP maps, so 
that £ Tl £ T2 = £ Tl -\- T2 . With this definition, Markovian 
quantum processes are analogous to their classical coun- 
terparts (see, for example, @, 0|), the classical require- 
ment of positivity being now replaced by complete posi- 
tivity as a result of the possible presence of genuine quan- 
tum correlations (entanglement) with some extra system. 
The structure of this kind of semigroups was analyzed in 
detail some time ago [§] , concluding that a quantum sys- 
tem will undergo a Markovian dynamics provided that 
its evolution satisfies a Master equation of the standard 



(Lindblad) form: 

| = £(p) = +£ 7k (vupVl - \{vlv kl p}) . 

k ^ ' 

(1) 

Here H is a self-adjoint operator, 7^ > 0, Vfe, and £ 
is called the generator of the semigroup (£ T — e £r ) . In 
analogy to the case of time dependent Hamiltonians in 
closed systems, the generators can be also time depen- 
dent Ct] then it is possible to prove that if and only if 
Ct can be written in the standard form ([1]), with H(t), 
7k{t) > and Vk(t) potentially time dependent, the fam- 
ily of propagators £(t 2 .t!) = Texp ( £ T drj satisfying 



2 



the composition law 



£(*2,t ) — %2,«i)£(ti,to)> 



(2) 



for all t2 > t\ > to > are CP maps. This situation is 
sometimes referred to as time-inhomogeneous Markovian 
dynamics. Note that the condition ([2]) is the quantum 
counterpart to the classic Chapman-Kolgomorov equa- 
tion. 

Conceptually, one could think of introducing a Marko- 
vianity measure via some type of optimization problem, 
such as evaluating a quantity of the form 



maxmin||£ (to+Mo) - £ M 



e>0 



(t +£,to)"' 



where || • || denotes some appropriate operator norm. The 
minimum is taken over the set of Markovian maps £ M , 
and the maximum over final times deals with the time 
continuous dependence of the dynamical maps (note that 
£(t ,t a ) — 1 independently of £ , and 1 is trivially Marko- 
vian). However, this quantity is hard to compute in prac- 
tice due to the nonconvex structure of the set of Marko- 
vian maps £ M 

In this Letter, we adopt a novel strategy and propose 
two possible ways to quantify the non-Markovian char- 
acter of a quantum evolution which avoid the definition 
of an optimization problem. Our key element will be ex- 
ploiting the specific behaviour of quantum correlations 
when a part of a composite system is subject to a lo- 
cal interaction that can be modeled as a trace-preserving 
CP map. This will allow us to give a necessary condi- 
tion to measure deviations from Markovianity even when 
the actual form of the dynamics is completely unknown, 
although some non-Markovian evolutions may be unde- 
tected. A necessary and sufficient condition can be pro- 
vided in the case when the dynamics is amenable to com- 
plete characterization, for instance, via quantum process 
tomography. Then, we are able to quantify strictly the 
Markovian character of the evolution and deviations from 
strict Markovianity can be unambiguously characterized. 

Let us consider first the case where we do not have 
any information about the dynamics of our system of 
interest, which we will consider initially to be a gen- 
eral, possibly composite d-dimcnsional quantum system. 
Our aim is to introduce a measure, that we denote by 
1^ , that quantifies the deviation from Markovianity in 
the evolution of the system. For that, we will initially 
prepare a maximally entangled state with an ancillary 
system which has to remain isolated from the decoher- 
ence sources, as illustrated in Fig. 1. Since local trace- 
preserving CP maps do not increase the amount of en- 
tanglement Q, it is evident from the composition law 
([2]) that the decay of the entanglement with an ancil- 
lary system will be monotonically decreasing for Marko- 
vian evolutions. This fact also prevents the formation 
of loops in diagrams concurence vs purity as illustrated 
in [10(. However, if the evolution is non-Markovian, 
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Parameters: 

-W =21, Af = 0.0209, 

- u= 10, 

- uj c = 15, w niax — 50, 

- Oscillators in Bath: 350, 

- Squeezing parameter r — 4. 
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FIG. 2: Results of the simulation for the non-Markovianity as 
a function of the strength of the coupling a between a damped 
harmonic oscillator and a bath with Ohmic spectral density 
J(cj) = aoje - "/" , for different temperatures. 



the requirement of strict monotonicity does no longer 



hold 
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as environmental correlations can lead to bi- 
partite entanglement to be increased and decreased as 
a function of time (as exemplified by the red curve in 
Fig. 1). Hence a conceptually simple way to quantify 
the degree of non-Markovianity of an unknown quantum 
evolution would be to compute the amount of entan- 
glement between system and ancilla at different times 
within a selected interval [iojimax] & n d check for strict 
monotonic decrease of the quantum correlations. That is, 
for AE = E{psA(t )} - E[p S A(t m ^)} (where E denotes 
some entanglement measure) and some initial maximally 
entangled system-ancilla state, |$) 
Pse{0) = | $)($|, we have 



1 T-^d-l 



Vd 



x{ E) 



r 


dE[ps A (t)} 


)t 


dt 



dt - AE, 



in such a way that if the evolution of the system is Marko- 
vian the derivative of E[psA(t)] is always negative and 
= 0. 

Note that since the knowledge of the exact form of the 
dynamics is not necessary to measure this method 

can be particularly useful in the study of infinite dimen- 
sional systems, where the computation of the exact dy- 
namical map is often difficult. For the sake of illustration, 
let us consider a single damped harmonic oscillator, with 
total system-bath Hamiltonian given by 



M 



M 



j'=i i=l 

where M is the number of oscillators in the bath, and we 
have assumed the validity of the rotating wave approx- 
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Parameters: 



- W =21, At -0.0417, 

- u = 10, 

- UJ C = 3, Wmax = 50, 

- Oscillators in Bath: 350, 

- Squeezing parameter r — 4. 
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FIG. 3: The analogous to Fig. 2 for a super-Ohmic spectral 



density J (to) = au> e 



note the different values of the 



cut-off frequencies and the order of magnitude of a. 



imation (RWA). The bath is assumed to be initially in 
a thermal state ps = exp (— Hb/T) /tr [exp (— Hb/T)], 

Hb = J2jLi UJ j a ] a j7 an d the system oscillator will 
be initially entangled with another oscillator, the an- 
cilla, in a two- mode vacuum squeezed state |A) = 
Vl — A 2 ^ n \ n )\ n )i where A = tanhr and r is the so- 

called squeezing parameter. Recall that this state is the 
most entangled Gaussian state at fixed mean energy n = 
sinh 2 r, and for infinite squeezing r — > oo it approaches to 
the maximally entangled state. So it seems appropriate 
to apply our method, although the shape of the curves 
does not depend significatively on r. Since the Hamil- 
tonian is quadratic in the annihilation and creation op- 
erators, it preserves the Gaussian character of the states 
and the amount of entanglement between system and an- 
cilla can be computed easily by means of the logarithmic 
negativity E N (p AB ) = log 2 ||p5s||i, where T A denotes 
the partial transposition with respect to the subsystem 
A and || • || i the trace norm 12j. To visualize the sensitiv- 



ity of the proposed measure I^ E \ two different spectral 
densities of the bath have been considered, as well as sev- 
eral initial temperatures. In Fig. 2, the behaviour of 
has been plotted for an Ohmic spectral density with ex- 



ponential cut-off J(w) = Ylj=i 9j — w i) 



-ul/uj r _ 
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where uj c is the so-called cut-off frequency. In order to 
make the analysis simpler, we focus only on the dynam- 
ical behaviour under different strengths of the coupling 
between system and bath. In the simulation, of course, 
the discrete and finite number of oscillators in the bath 
will affect the results, but in order to avoid this effect we 
have taken a number of oscillators distributed in a range 
of frequencies (see details on the plots) such that our fi- 
nal time i max is shorter than the recurrence time of the 
bath but large enough to capture the general properties 



of the dynamics. It is clear from Fig. 2 that in the limit 
a — > we approach a Markovian evolution independently 
of the value of T. This is the well-known weak-coupling 
limit, whereas the on-set of deviations from Markovian- 
ity for stronger coupling strengths are T-dependent, so 
that the value of a for which the measure l( E > becomes 
nonzero increases with increasing temperature, the same 
happens with the value of I^ E ' itself for a large enough. 
An analogous situation is encountered in the presence 
of a super-Ohmic spectral density J {to) = au'e""''" 1 , 
as shown in figure 3; however, the specific value of 1^ 
for given T and a depends strongly on the bath spectral 
function. 

There can be however non-Markovian quantum evolu- 
tions that remain undetected by the proposed measure, 
given that the requirement X^> > is a sufficient con- 
dition for deviation from Markovianity. A necessary and 
sufficient condition can nevertheless be formulated if the 
specific form of the quantum evolution, as given by some 



dynamical map £ , 



(Mo) 



is amenable to exact reconstruc- 



tion between the some initial time to and a final time t. 
This can in principle be done by means of process tomog- 
raphy or perhaps resorting to a theoretical microscopic 
model. Let us take our initial time as t — without loss 
of generality. Then, because of the continuity of time, we 
can split the dynamical map as 



£(*+£, o) — <%+£, t)£(*,o) 



(3) 



for any times t and e. If the time evolution implemented 
by £(t,o) is Markovian, we have already mentioned that 
£{t 2 ,tx) is CP for any intermediate times (t + e) > f 2 > 
t\ > (that is, £(t,o) is infinitesimally divisible in the 
sense of [j| for any t). However if and only if there exist 
times t and e such that £(t+ e ,t) is not CP, the dynamics 
will be non-Markovian. Note that this is the ultimate 
reason behind the possible increase of the system-ancilla 
entanglement at some local times, which is the basis of 
our previous measure. This partition can be extracted 
from the known dynamical map £(t,o) as a function of 
t just by applying £(7o) ( wn i cn may not be a CP map) 



to the Eq. ((3]). Since £(t^>o.o) —> 1, f° r t small enough 
£tt,o) w i u b e invertible, if for larger times £(t,o) is not 
invertible, we do not have enough information to define 
£(t+c,t) i n an unequivocal way (this is one consequence of 
being blind to one part of the whole system) . Then there 
are several routes to follow depending on the nature of 
the singularities of £(t,o)> f° r example, strategies based 
on pseudoinverse maps have been recently applied in a 
similar context [l3j |. 

Being |$) a maximally entangled state of our open 
system and some ancillary one, because of the Choi- 
Jamiolkowski isomorphism 1J| £^ t+e t ^ is CP if and only 
if (£( t +e,t) ® 1) |$)($| > 0. Hence, given the trace- 
preserving property, we can take the following definition 



4 



as a measure of the non-CP character of £r t+e>t ), 

f NCP (t + e,t)= || (£ (t+M) ® 1) (|$)($|) || x . 

Therefore, £ (t+e.t) is CP if and only if fNCp{t + e ,t) = U 
otherwise /jvcp(* + e, t) > 1- Now fNCp(t + £,t) will be 
the building block of our measure of non-Markovianity. 
To construct it we leave e to be infinitesimal to define the 
(right) derivative of fNCp{t + e , t)'- 



git) 



lim 



f NC p(t + e,t)-l 



noticing that g(t) > 0, with g(t) — if and only if £ (t+e,t) 
is CP. Therefore the integral 



g(t)dt 



can be taken as a measure of non-Markovianity, and as 
long as g{t) decreases fast enough (this will not be always 
the case) will be finite. Actually a normalized version of 
this measure can be T>nm = xTT> m sucn a way that 
"Dnm = for 1=0 (i.e. Markovian evolution) and 
T^nm — > 1 for 1 — > oo . 

To illustrate the behaviour of I. let us consider the 
dynamics of one qubit modeled by a possibly non- 
Markovian differential master equation = C t ip) 0, 
since in the limit e — > the solution of this equation 

aCt<L 15[, we only need to 



formally tends to £( t+e t ) — > e" 
expand it inside of trace norm up to first order to calcu- 
late git), 



[l + (£ t ®l)e]|$)($| 



1 



git) = lim 

For instance, for a very simplified evolution of a qubit 
such as pure dephasing written like & = j{t)io z po z — p) 
we immediately obtain 



<?(*) = 



and finally 



for 7(t) > 

-2j(t) for j{t) < 



-2 / j(t)dt. 

h(t)<o 



Therefore, for this kind of evolution, X is proportional to 
the area of jit) which is below zero. In particular, if the 
negative values of jit) do not tend to zero fast enough 
I — > oo [This is actually the case in the second example 
of d where j(t) ~ tan(t)] 0. 

In summary, we have proposed two different 
approaches for the problem of quantifying non- 
Markovianity of general quantum evolutions; one is based 
on a sufficient condition whose evaluation does not re- 
quire any prior knowledge of the quantum evolution itself, 
and the other provides a measure which quantify strictly 



the non-Markovianity, provided that the structural form 
of the dynamical map is known. Remarkably, the evalua- 
tion of the proposed measures does not require solving an 
optimization problem and would be suitable for provid- 
ing information on the deviations from Markovianity in 
the experimental implementation of effective spin models 
using controllable systems, as for instance, trapped ions 
17]. This type of experiments can prove extremely valu- 
able in the subsequent formulation of detailed models of 
system-environment coupling in complex systems, both 
in condensed matter and, potentially, in some biological 
aggregates. These compounds have recently become ex- 
perimentally probable at the femtosecond scale [l8| but, 
given their complexity, no detailed microscopic models 
for the interaction with their surroundings are currently 
available. 
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